arXiv:1508.06283vl [astro-ph.GA] 25 Aug 2015 


Draft version August 27, 2015 

Preprint typeset using DTeX style emulateapj v. 08/22/09 


STAR-FORMING BRIGHTEST GLUSTER GALAXIES AT 0.25 < ^ < 1.25: A TRANSITIONING EUEL SUPPLY 

M. McDonald^ B. Stalder^, M. Bayliss^’^, S. W. Allen^’®’^, D. E. Applegate®, M. L. N. Ashby®, M. Bautz®, 

B. A. Benson^’®®’®®, L. E. Bleem®®’®^’®®, M. Brodwin®^, J. E. Carlstrom®®’®®’®^’®®’®®, I. Cmu®®, S. DesaP®’®^, 
a. H. Gonzalez®®, J. Hlavacek-Larrondo®^, W. L. Holzapfel^®, D. P. Marrone^®, E. D. Miller®, 

C. L. Reichardt^^, B. R. Saliwanchik^®, A. Saro®®, T. Schrabback®, S. A. Stanford^^ A. A. Stark®, J. D. Vieira^®, 

a. Zenteno^® 

Draft version August 27, 2015 

ABSTRACT 

We present a multi-wavelength study of 90 brightest cluster galaxies (BGGs) in a sample of galaxy 
clusters selected via the Sunyaev Zebdovich effect by the South Pole Telescope, utilizing data from 
various ground- and space-based facilities. We infer the star formation rate (SER) for the BGG in 
each cluster, based on the UV and IR continuum luminosity, as well as the [O ll]AA3726,3729 emission 
line luminosity in cases where spectroscopy is available, finding 7 systems with SER > 100 M© yr“®. 

We find that the BGG SER exceeds 10 Mq yr“® in 31 of 90 (34%) cases at 0.25 < 2 : < 1.25, compared 
to ^1-5% at z ^ 0 from the literature. At 2 : > 1, this fraction increases to 92^.11%^ implying a steady 
decrease in the BGG SER over the past ^9 Gyr. At low-^, we find that the specific star formation 
rate in BGGs is declining more slowly with time than for field or cluster galaxies, most likely due to 
the replenishing fuel from the cooling IGM in relaxed, cool core clusters. At 2 : > 0.6, the correlation 
between cluster central entropy and BGG star formation - which is well established at ^ ^ 0 - is 
not present. Instead, we find that the most star-forming BGGs at high- 2 : are found in the cores of 
dynamically unrelaxed clusters. We investigate the rest-frame near-UV morphology of a subsample 
of the most star-forming BGGs using data from the Hubble Space Telescope^ finding complex, highly 
asymmetric UV morphologies on scales as large as ^50-60 kpc. The high fraction of star-forming 
BGGs hosted in unrelaxed, non-cool core clusters at early times suggests that the dominant mode 
of fueling star formation in BGGs may have recently transitioned from galaxy-galaxy interactions to 
IGM cooling. 

Subject headings: galaxies: clusters: general - galaxies: clusters: intracluster medium - galaxies: 

elliptical and lenticular, cD - galaxies: starburst - X-rays: galaxies: clusters 
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1. INTRODUGTION 

One of the great mysteries in astronomy today is why 
>90% of the baryons in the Universe, which a re in diffuse 
gas w ith relatively short cooling times fe.g.. iShull et al.l 
SqU, have not cooled and formed stars. This inefficient 
star formation manifests as a significant disagreement 
between the predicted galaxy luminosity function from 
AGDM cosmological simulations and that observed in 
the local Universe. In the most massive galaxies in the 
Universe, which are found at the centers of rich galaxy 
clusters, this disagreement is maximized, with central 
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cluster galaxies being substantially less m assive than pre¬ 
dicte d by simple models (see review by iSilk fc Mamonl 
l2Q12f ) . This has become known as the “cooling flow prob¬ 
lem” and can be stated simply as: “Why, given the short 
cooling time of the intracluster medium in the cores of 
some galaxy clusters, do we not observe massive star- 
burst galaxies at the centers of these clusters?”. 

Over the past couple of decades, much effort has 
been devoted to answering this question. Early 
work focused on searching for multiphase gas and 
star formation in brightest cluster galaxies (BCGs). 
Numerous studies have found evidence for ultra- 
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tive of ongoing or recent star formation. Star-forming 
BCGs were found preferentially in galaxy clusters 
with “cool cores”, as identified by a central density 
enhancement in t h e ICM (e.g., IVikhlinin et al.l I2QQ7I: 
iSantos et al.l I2QQ8I: [Hudson et al.l 1201011 or low central 
entropy/cooling tim e fe.g., iCavagnolo et al.ll2QQ8L I2QQ9I: 
[Hudson et al.ll^CUOl ). These and other works established 
a link between the cooling ICM and the presence of 
multiphase gas, suggesting that cooling flows may 
indeed be fueling star formation in BCGs. However, 
the typical star formation rates inferred from a vari¬ 
ety of indicators were found to b e only ^1% of th e 
expected ICM cooling rate (e.g., [Q’Dea et al.l [2QQ8f ). 
Roughly one of the two orders of magnitude in this 
disagreeme nt can be accounted fo r by inefficient star 
formation (jMcDonald et al.l l2Q14aI ) , however a further 
order-of-magnitude disagreement between the cooling 
predictions and observations still remains. 

Some form of feedback is necessary to prevent the bulk 
(^90%) of the cooling ICM from becoming fuel for star 
formation. The leadi ng candidate is “radio-mode” feed¬ 
back (see reviews bv lFabianll2Q12l: iMcNamara fc NulsenI 
I2Q12I ) from active galactic nuclei (AGN), which are ubiq¬ 
uitous at the centers of cool core clusters (|Sunll2QQ9l) . The 
mechanical energy output from these AGN are sufficient 
to offset cooling on large scales, preve nting runaway cool¬ 
ing in the maj o rity o f clusters (e.g., iBirzan et al.ll2QQ^ 


rio rity o i clusters (e.g., irza 
Rafferty et alT I2QQ8I: iHlavacek-Larrondo et al.l 12012 
20141 ) with a few notable exc eptions ([McNamara et al" 


2QQ6I : [McDonald et al.l l2Q12c[ ) . The low levels of star 


formation and gas in multiphase filaments are under- 
stood to be local thermodynamic in s tabilities (e.g., 
Sharma et al.ll2QlQl: iCaspari et al.ll2Q12l: iMcCourt et al.l 


20121: I Volt fc Donahuel I2015D in regions where, locally^ 


cooling dominates over feedback, despite the overall 
global balance. These star formation rates, which av¬ 
erage a few Mq yr“^, may contribute a few percent to 
the total stellar mass of the BCG over the past ^8 Gyr - 
the majority of the growth in these systems likely comes 


from “dry mergers” (mergers of gas-poor galaxies), which 
increa se the stellar mass by a factor of ^2 from z = 1 to 
z = 0 (|Ruszkowski fc Springelll^OOOD . 

The picture presented here is based almost entirely 
on observations of nearby {z < 0.3) galaxy clusters. 
This is, in part, due to the fact that these clusters are 
more easily studied because of their proximity (improved 
signal-to-noise, angular resolution, etc). Equally impor¬ 
tant, however, is the scarcity of well-understood sam¬ 
ples of high-redshift galaxy clusters. Until recent years 
there were few samples of galaxy clusters with known 
masses at 2 : ^ 0.5 - surveys that did probe high red- 
shift clusters were generally flux-limited or were assem¬ 
bled from serendipitous detections via a heterogeneous 
collection of methods. With the advent of large-area mm- 
wave surveys utilizing the Su nyaev Zel’dovich (SZ) effect 
([Sunvaev fc Zeldovichl 119721 ) to detect galaxy clusters, 
this situation has changed dramatically over the past 
several years, with the latest surveys achieving nearly 
redshift-independent detection of clusters above a fixed 
mass threshold at 2 : > 0.3. Most recently, the com¬ 
pleted 2500 deg^ SZ survey with the South Pole Tele¬ 
scope fSPT: ICarlstrom et al.|[2QTTl) has discovered more 
than 500 mass ive galaxy clusters , the majority of which 
are at 2: > 0.5 (|Bleem et al.ll201§ ). In this work, we focus 
on a subsample of this survey which has b een targeted 
for X-ray follow-up (|McDonald et al.|[2013bD and spans a 
redshift range of 0.3 < 2: < 1.2. The availability of opti¬ 
cal photometry and spectroscopy for the majority of the 
BCGs in this sample, along with archival UV (GALEX) 
and IR (WISE) data allows us to study star formation 
in BCGs at high redshift for a complete, mass-selected 
sample. 

The remainder of this paper is structured as follows. In 
§2 we define the sample used in this work, and present 
the multiwavelength data and analysis techniques that 
will enable us to identify star-forming BCGs. In §3 we 
isolate the sample of star-forming BCGs and attempt to 
determine whether there is any evolution in their proper¬ 
ties or the properties of their host clusters. In §4 we dis¬ 
cuss these results, addressing bias and selection concerns, 
while trying to draw a broad picture of galaxy and galaxy 
cluster evolution within which these results fit. We finish 
in §5 with a brief summary of the important results of 
this work, and a look toward the future. Throughout this 
work we assume Hn = 70 km s~^ Mpc“^, Um = 0.27, 
Oa = 0.73, and a iSaloet^ (|1955l ) initial mass function 
(IMF). 


2. DATA & ANALYSIS 
2.1. Cluster Sample and BCG Seleetion 

W e initially define the samp le to include all 83 clusters 
from lMcDonald et al.l (|2013bl ). which were selected from 
the SPT 2500 deg^ survey (|Bleeni et al.ll201^ and sub¬ 
sequently observed in the X-ray with the Chandra X-ray 
Observato ry. To this sample, we add an additional 8 clus¬ 
ters from iBleem et al.l (|2015f ) that have archival Chan¬ 
dra data from other sources (SPT-CLJ0106-5943, SPT- 
CLJ0232-4421, SPT-CLJ0235-5121, SPT-CLJ0516-5430, 
SPT-CLJ0522-4818, SPT-CLJ0658-5556, SPT-CLJ2011- 
5725, SPT-CLJ2332-5053). X-ray data products, such as 
central entropy (Kp), are presen t ed for this full sample 
of 91 clusters in lMcDonald et all (|2013b[ ) - we direct the 
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reader there for a full description of our X-ray method¬ 
ology. 

For each cluster, we have obtained some combination 
of ground-based g^r^i^z optical imaging, ground-based J, 
H, K near-IR imaging, and Spitzer IRAQ 3.6/im , A.bpm 
imaging, as described in detail in iBleem et al.l (j2Q15f ). 
These data have been used to provide optical confirma¬ 
tion of the SZ-selected clusters, and to estimate a photo¬ 
metric redshift based on the member galaxy colors. We 
remove from our sample all clusters for which the optical- 
IR follow-up was performed in <2 filters, preventing a 
reliable fit to the stellar continuum. After this cut, we 
were left with a sample of 82 clusters. From this follow¬ 
up imaging, we select the BCG^ initially as the brightest 
red-sequence galaxy within an apertur e of Ro od centered 
on the SZ peak, following iSong et all ^2012\ ). We then 
visually inspect each cluster and select a new BCG if 
one, or both, of two conditions are met: i) there is a 
similarly-bright galaxy that is significantly closer to the 
X-ray peak (29% of systems); ii) there is a bright blue 
galaxy on or near the X-ray peak that was not selected 
due to our initial preference for red galaxies (3% of sys¬ 
tems). Wherever possible, we use HST imaging from the 
SPT weak-leasing follow-up programs with IDs 12246, 
12477, 13412 (Pis: Stubbs, High, Schrabback) to aid in 
the visual identification of the BCG. 

In all cases where it was not clear which galaxy was 
the BCG, we retain the original red BCG, allowing for 
two BCGs in these clusters. In this case, we give each 
BCG a 50% weight when performing statistical analysis 
and include the limiting cases in our error estimates (i.e., 
upper limit contains star-forming BCG, lower limit con¬ 
tains passive BCG). This procedure resulted in a sample 
of 90 BCGs in 82 clusters. The net effect of allowing mul¬ 
tiple BCGs in cases where identification of a single object 
was challenging is to increase our uncertainties, increas¬ 
ing the likelihood that the “true” answer lies without our 
error bars. 

2.2. X-ray Analysis: Central Entropy and Luminosity 

Several studies have found correlations between the 
amount of star formation in the BCG and the 
core e ntropy (e.g . . iCayagnolo et al.l 120081: 10 ’Pea et al.l 
120081: IVoit et alT 12015) and luminosity-derived cool¬ 
ing rate fe.g.. ICrawford et al.l 1 19991 : 10 ’Pea et al.l 120081: 
iMcDonald et al.l l2010L l2011bf ). This has led to the 
conclusion that star formation i s being fueled i n the 
BCG by residual cooling flows (IVoigt fc FabianI l2004 
iTremblav et al.l 120121: iMcDonald et al.l l2014a[ ). In or¬ 
der to test whether this trend was established at 
high-^, we require estimates of the core entropy and 
cooling rate for each cluster. Given that we only 
have ^2000 X-ray co unts per cluster , mod eling the 
central entropy (e.g., iCavagnolo et al.l l2009[) or esti- 
mating the spectrosc opically-derived cooling rate (e.g., 
IVoigt fc FabianI [20041) is not feasible. Instead, we com¬ 
pute spectroscopic quantities (bolometric luminosity, 
temperature) from a circular aperture with radius of 

^ We note here that the term “BCG” is largely inappropriate for 
this work, since there may be brighter galaxies within the virial ra¬ 
dius (e.g., AGN). The more appropriate monicker would be “central 
cluster galaxy”, but we opt for BCG throughout this paper since 
it is more commonly used. 


0.075 R 500 (where R 5 nn was derived ba sed on the Yx- 
M 500 relation of IVikhlinin et al.l l2009[ ) , which should 
roughly corres pond to the deprojecte d core tempera¬ 
ture (see e.g.. IMcDonald et al.M2014b[) . X-ray spectra 
extracted from this aperture are modeled with a pho¬ 
tometric absorption (PHABS) and plasma (apec) model, 
allowing the temperature, metallicity, and normalization 
of the plasma model to vary. This choice of aperture 
is meant to capture the “core” properties, reflecting the 
realistic fuel reservoir that the BGG may have access to. 

The “central” density is computed based on the 
de projected X-ray surfa ce b rightness profile, follow¬ 
ing IVikhlinin et al.l ()2QQ6f ) and IMcDonald et al.l (j2Q13b[ ) . 
Since the measurement of electron density requires far 
fewer X-ray counts than the measurement of spectro¬ 
scopic temperature, the central density is measured at 
r = O.OIR 500 , or roughly 10 kpc for a typical cluster in 
this sample. We combine the projected core tempera¬ 
ture and deprojected core density to arrive at a pseudo- 

deprojected core entropy (Kq = kTo,or 5 R^^^ 

Unlike our previous work (jMcDonald et al.l l2Q13b[ ). 
here we measure “central” quantities at the X-ray peak, 
rather than the large-scale centroid. In general, the 
BCG is located closer to the X-ray peak than the cen¬ 
troid in cases of merging clusters, which motivated this 
choice. We will discuss in later sections the effects of this 
choice. For a more detailed description of our X-ray anal¬ 
ysis tec hniq ues, we direct the reader to IMcDonald et al.l 
(|2Q13b[ ) and IMcDonald et all (j2Q14b[ ). 

2.3. UV-Optical-IR Photometry and SED Modeling 

Ground-based optical and/or near-infrared imaging for 
all clusters in this sample have been obtained at a variety 
of wavelengths as part of a confirmation and photomet¬ 
ric redshift follow-up campaign. The acquisition, reduc¬ 
tion, and calibrat i on of thes e data are pr e sente d in de¬ 
tail in iSong et al.l (j2Q12[ ) and iBleem et al.l (j2Q15[ ). Aper¬ 
ture pho tometry for the BCG w as obtained from SEx- 
tracto r (iBertin fc ArnontsI Il996f) . following ISong et all 
(|2Q12f ) and iBleem et al.l (j2Q15f ) . 

In addition to these existing data, we have ac- 
qnired new ground-b ased n-band imaging with Megacam 
(jMcLeod et al.l[2QT^ on the Magellan Clay telescope for 
49 clusters in this sample. These 49 clusters were those 
that had the least restrictive upper limits on the BCG 
SFR at other wavelengths, generally lacking in spectro¬ 
scopic or deep GALEX coverage. The exposure time for 
these observations was chosen to provide an overall sen¬ 
sitivity of our survey to obscured star formation rates 
of ^10 Mq yr“^ - without these additional data, our 
sensitivity limit would vary significantly with redshift. 
These data were reduce d using the s t andar d photomet¬ 
ric pipeline described in IBleem et al ] (l2nT5h . 

The position of each BCG was cr oss-referenced 
with the GALEX‘S (|Morrissev et al.l 120071 ) and WISE^ 
(| Wright et al.ll2010[ ) archives, assuming a maximum off¬ 
set of 2", from which we obtained near-UV (NUV) and 
near-mid IR (NIR, MIR) photometry for each BCG. If 
the BCG was undetected by either of these surveys, we 
instead obtained upper limits. Spitzer 3.6/im, 4.5/im, 

^ http://galex.stsci.edu/GR6/ 

^ http://irsa.ipac.caltech.edu/Missions/wise.html 























































































4 


and WISE 4-band photometr y was conver t ed fro m Vega 
to AB magnitudes following iPrice et ahl (j2QQ4f ) . The 
resulting UV-optical-IR spectral energy distributions 
(SEDs), which span 2000A to 22 /im in the observed 
frame, are shown in Figure lA.ll for our sample of 90 
BCGs. As mentioned above, we have discarded all BCGs 
for which we have <2 photometric measurements (ex¬ 
cluding upper limits). 

The observed SEDs are fit in two stages. First, we 
model the full SED with a single-age population with 
formation redsh ift z/, solar metallicity, and a Salpeter 
(jSalpeterl Il955[ ) initial mass function. We assume a 
uniformly-distributed range of 2 ;/ from 2-5 in our mod¬ 
els, which leads to some uncertainty in our resulting stel¬ 
lar masses and star formation rates. The spectrum for 
this old s tellar population was generated using Star- 
BURST99 (jLeitherer et al.l[l999l ) and was convolved with 
our broadband filter set. The model spectrum was fit 
to the data using mpfitfun^, which minimizes with 
respect to the two free parameters (normalization, red- 
shift). The redshift was allowed to vyy with in the mea¬ 
sured uncertainty from iBleem et al.l ()2Q15f ). while the 
normalization, which corresponds to the stellar mass of 
the old population, was left free. 

Beyond this single-component model, we also consider 
additional components in the UV and IR representing 
contributions from a young stellar population and warm 
dust, respectively. At short wavelengths, we model ex¬ 
cess emission using Starburst99, assuming a constant 
star formation rate over the past 30Myr (roughly the 
AGN duty cycle). We note that adjusting this timescale 
down to lOMyr or up to 100 Myr results in deviations in 
the derived SFR of ^20%. We assume that the emission 
from young stars is obscured by dust, incorporating the 
iCalzetti et al.l (|2000[) extinction law and an intrinsic red¬ 
dening of E{B—V) = 0.3T0.1 - this rang e is based on ob¬ 
serva t ions of nearby star-form ing BGGs (|Grawford et al.l 
1 19991: [McDonald et al.l l2012aD . The young stellar com¬ 
ponent has a single free parameter, corresponding to the 
ongoing, extinction-corrected star formation rate. At 
long wavelengths, we mimic a dusty component with a 
mid-in frared power-law {Fy oc following iGasevI 

(|2Q12[) . This dust component is artificially truncated at 
shorter wavelengths, so that it will not add UV flux. We 
only include the young and dusty components if their 
inclusion improves the which is the case for 34% 

(31/90) of the BCGs in our sample. For the remaining 
66% (59/90), the “passive-evolution” model yields a suit¬ 
able fit to the d ata. The results of this SED-fitting are 
shown in Figure IAtI 

2.4. Optical Spectroscopy 

For 36 of the 95 BCGs in this sample, we have op- 
tical spectroscopy f rom a combination of the IMACS 
([Dressier et al.|[2Mll) and LDSS3 spectrographs on Mag¬ 
ellan. These spectra were initially obtained as part of 
a spectroscopic redshift campaign which targeted, pri¬ 
marily, red sequence galaxies. The full details of this 
spe ctroscopic fol l ow-up is provided in iRnel et ahl (|2Q13D 
and IBleem et al.l (|2Q15D . For each spectrum, we measure 
the [O ll] and 116 equivalent widths and the 4000A break 

http://www.physics.wisc.edu/~craigm/idl/down/mpfitfun.pro 


streng th (D 4 ooo )5 via indices defined by iBalogh et al.l 
(IT^ . [On] equivalent widths are converted to emis¬ 
sion line fluxes using the continuum level, as determined 
by the SED modeling described in §2.2, at rest-frame 
3727 A. From these 36 spectra we find significant (>3cr) 
[O ll] e mission in 5 systems, one of which is the Phoenix 
cluster (|McDonald et al.ll2Q12cD . This system allows us 
an opportunity to compare the index-based measurement 
of [On] t o the total, flux from de ep integral-field spec¬ 
troscopy (|McDonald et al.l l2Q14cl) . Using the long-slit 
data, the index-based technique of IBalogh et al.l (|1999f ) 
combined with our best-fit optical SED yields a line flux 
of f[oii] = 2.2 X 10 ergs s ^ cm ^ for the central 
galaxy in the Phoenix clu s ter, wh ile the total measured 
flux from [McDonald et al.l (|2Q13a[) is f[oii] = 1.6 x 10“^"^ 
ergs s“^ cm“^. Thus, for this single comparison, the 
index-based approach does an adequate job of reproduc¬ 
ing the measurement made via emission-line modeling of 
significantly higher-quality data. 

When correcting the [O ll] emission line flux for extinc¬ 
tion, we assume that the ionized gas has a factor of two 
higher redden i ng th a n the ste l lar co ntinuum, following 
iGalzetti et al.l (|2QQQD : iGalzettil (|2QQlD . This assumption 
is motivated by the idea that the warm ionized gas traces 
the highly-ionizing O and B stars, which tend to be em¬ 
bedded in dusty, star-forming regions. 

2.5. Star Formation Rates and Their Uncertainties 

2.5.1. UV-Derived SFRs 

The calculation of star formation rates (SFRs) from 
the UV photometry is straightforward. In §2.3 we 
describe our SED-fitting procedure, which includes an 
intrinsically-reddened, continuously-star forming popula¬ 
tion derived from Starburst99 models. The normaliza¬ 
tion of this component yields the current star formation 
rate, under the assumption of constant star formation for 
the past 30 Myr. 

The uncertainty in the UV-derived SFR is dominated 
by our uncertainty in the amount of intrinsic extinction 
{E{B — V)) and, to a lesser extent, the formation redshift 
of the old stellar population (zf). In order to propagate 
our uncertainty in these quantities to our estimate of the 
SFR we perform 100 fits to each SED, varying Zf and 
E{B — V) in each fit. We assume a normal distribution 
for E{B—V)^ with E{B—V) = 0.3T0.1, i notivated by ob¬ 
serva t ions of nearby star-for ru ing BGGs ([Grawford et al.l 
1 19991: [McDonald et al.|[2012aD . For the old stellar pop¬ 
ulation, we assume a uniform distribution of formation 
epochs, from 2 :/ = 2.0 — 5.0 Gyr. This Monte-Carlo ap¬ 
proach results in typical uncertainties in the UV-derived 
SFRs of 0.38 dex, or a factor of ^2.5. 

2.5.2. [Oil]-Derived SFRs 

As described in §2.3, [Oil] fluxes ye de rived from a 
combination of indices (jBalogh et al.lll99^ . which pro¬ 
vide an estimate of the equivalent line width, and SED 
fitting, which provides the interpolated continuum level 
at 3727A. We convert the measured emission line flux 
to a SFR, assuming SFR[o 77 ] = 9.53 x 10“^^ L[o 77 ] 
(jKewlev et al.l[2003 ). As in §2.5.1, we assume a normal 
distribution of intrinsic r eddening, with E(B — V)nas = 
2E{B — V) stars, following ICalzetti et al.l (|2000f ): IcSzettil 
(l200lh . As is the case for the UV-derived SFRs, our 
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uncertainty in this extinction corrections dominates the 
uncertainty in the inferred SFR. Propagating this uncer¬ 
tainty through the fitting procedure via a Monte Carlo 
approach yields typical uncertainties of 0.19 dex, or a 
factor of ^1.5. 


2.5.3. IR-Derived SFRs 

In the mid-infrared, we derive the SFR by first extrap¬ 
olating our power-law fit from observed 22jam to rest- 
frame 24/im. The uncertainty on this calculation is both 
redshift dependent, since at higher redshift we are ex¬ 
trapolating over a larger wavelength range, and depen¬ 
dent on the assumed power-law slope. We incorporate 
the uncerta inty in the power-law slope {a = 2.0 ± 0.5; 
lCasevll^012f ) in our calculation of the 24/im luminosity by 
assuming a normal distribution of values and performing 
100 fits to the data. The average resulting uncertainty in 
L 24 /im is 0.2 dex, or a factor of ^1.6. These extrapolated 
values of L 24 ^rn ar e then converted into an estimate of 
the SFR, following iCalzetti et al.l (j2007f ). 

2.5.4. Comparison of SFR Estimates 

In Figure [U we compare SFRs derived from the three 
indicators described above. We have only limited overlap 
between the three subsamples, with only 3 BCGs hav¬ 
ing both IR- and UV-derived SFRs and an additional 
4 BCGs with both UV- and [O ll]-derived SFRs. Only a 
single BCG (SPT-CLJ2344-4243) has SFRs inferred from 
all three methods. Nonetheless, we proceed to compare 
how well these estimates agree for systems in common. 
In general, SFR estimates from different indicators agree 
within the systematic errors. Assuming a one-to-one re¬ 
lation, we measure a scatter of 0.37 dex between UV- 
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and IR-derived SFRs, and 0.15 dex between UV- and 
[O ll]-derived SFRs. The better agreement between the 
UV- and [O ll]-derived SFRs is most likely due to the fact 
that the dominant systematic uncertainties (extinction, 
old stellar population age) are correlated between these 
SF R est imates. SFRs for each BGG are listed in Ta¬ 
ble lA.li along with relevant information about the host 
cluster. 


2.6. AGN Contamination 

Before presenting results from this survey, we would 
like to draw attention to one potential issue with our ap¬ 
proach. All three indicators of star formation employed 
here - UV continuum, [O ll] line emission, and IR contin¬ 
uum - are also indicators of active nuclei. The relative 
amounts of contamination at each of these wavelengths 
depends on the type of AGN considered. For example, 
in the Phoenix cluster, more than half of the total IR 
continuum comes from a dusty QSO. On the other hand, 
this AGN contributes <5% to the UV and [Oi l] emis¬ 
sion. Several of the SEDs shown in Figure lA.II exhibit 
a powerlaw shape that is consistent with both a dusty 
starburst and an AGN. 



-UO -0,5 0,0 0,5 

[3.5] - [4.5] (AB) 

Fig. 2.— Distribution of mid-IR color, [3.5] — [4.5], for the BCGs 
in this study from the Spitzer (red) and WISE (blue) space tele¬ 
scopes. The vertical cut at [3.5] —[4.5] = 0, wh ich separates galaxies 
from AGN, is motivated bv I Stern et al.l (| 20051 b and is valid for sys¬ 
tems at 2 : < 1.3. This figure demonstrates that, with the exception 
of the Phoenix cluster (SPT-CLJ2344-4243), this sample of BCGs 
is relatively free of strong AGN. 


10 100 1000 
SFR|r , SFR[o II] 

Fig. 1. — Comparison of SFRs derived via UV, [On], and IR re¬ 
lations for all 6 systems with multiple measurements. The dashed 
line represents the one-to-one relation, while the dotted lines show 
the measured scatter between the UV-IR (red) and UV-[0 ll] 
(blue) relations. The dominant source of uncertainty for UV- and 
[O ll]-derived SFRs is the extinction correction, while for IR it is 
the extrapolation from observed-frame 22pm. to rest-frame 24p,m. 
This figure demonstrates the overall agreement between these dif¬ 
ferent indicators, despite the significant uncertainties involved in 
estimating the SFR. 


In an attempt to quantify the contamination due to 
AGN in this sample, we consider in Figure [2] the 3.5/im - 
4.5/im color for 49 BCGs with Spitzer data a nd 82 BCGs 
with WISE data. Color corrections from iStern et ahl 
(12011 ) have been applied to match data from these tele- 
scopes to a common photometric system. Following 
IStern et a] J (l 2 fiM l 2 (m . we classify systems with mid-IR 
colors [3.5] — [4.5] > 0 as AGN, while those with redder 
colors are either passive or star- forming galaxies. In this 
color space, the Phoenix cluster (jMcDonald et al.ll2QI2c[ ) 
is the only BCG harboring a strong AGN. At most, we es¬ 
timate that ^4 clusters in this sample may harbor strong 
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AGN at their centers, based on this mid-IR color selec¬ 
tion. This lack of AGN contamination is further con¬ 
firmed by visual inspection of the Chandra data for each 
cluster, which show a general lack of strong X-ray point 
sources coincident with the BCG in all clusters with the 
exception of Phoenix. 

We proceed with this work assuming that all systems 
identified as star-forming are, indeed, star-forming, but 
remain cognizant of the fact that some fraction of these 
may host an AGN. We will return to this dilemma in the 
discussion section. 

3. RESULTS 

Below, we summarize the key results to emerge from 
this data set. We defer a detailed discussion of these 
results to the discussion section. 

3.1. Star Formation Rates in BCGs at 0.25 < z < 1.2 

In Figure [3] we show the star formation rates, derived 
via three different methods, for the sample of BCGs de¬ 
scribed in §2.1. This plot demonstrates that the WISE 
mid-IR data typically provide the best upper limits at 
low-redshift, but are unable to provide meaningful upper 
limits on the star formation rate for BCGs at 2 ; > 0.5. 
At 2 ; > 0.5, the most restrictive limit generally comes 
from our deep i4-band follow-up program, which was de¬ 
signed to achieve a redshift-independent sensitivity. For 
the most part, we are sensitive to SFRs higher than 10 
Mq yr“^, with <10% of systems having limits higher 
than this threshold. 

In the right panel of Figure [3l we combine the con¬ 
straints from the three different SF indicators. This 
plot shows a significant number of BCGs with SFRs 
from 10-300 Mq yr“^ at 2 : > 0.4. One quarter (7/31) 
of the star-forming BCGs have SFR >100 Mq yr“^ - 
of these, 5 are at 2 : > 0.9. For comparison, prior to 
this work there were only four confirmed clusters with 
extinction-co rrected SFR >100 Mr?) yr“^ in their BCG: 
Abell 1835 (iMcNamara et al.l [2003 l. RX J1504.1-0248 
(IQgrean et al.ll2010F MACS J1 931.8-2634 dEhlert et all 
l201l[ ). and the Phoenix cluster (jMcDonald et al.ll2012c[ ). 
The median star formation rate for the 31 star-forming 
BCGs identified over the full redshift range is 50 Mq 
yr“^. It is worth noting that 5 of the 7 most star-forming 
BCGs (SFR > 100 Mq yr“^) are at 2 : > 0.9, despite the 
fact that only 15% of the clusters in this sample are at 
such high redshift. Likewise, 3 of the 5 most star-forming 
BCGs are at 2 : > 1.1, despite this representing 6% of the 
sample. 

We also show in the right panel of Figure [3] the aver¬ 
age SFR as a function of redshift in five redshift bins. 
These averages, provided in Table [U have uncertainties 
that are derived by bootstrapping errors on individual 
data points. For non-detections, we assume zero SFR for 
the lower-bound of the uncertainty and the upper limit 
for the upper bound. Thus, the ranges shown encom¬ 
pass both the statistical uncertainty in the mean from 
detections and the uncertainty in the true value of the 
non-detections. This analysis reveals a strong evolution 
in the average star formation rate in BCGs over the past 
9 Gyr. This increase in the average SFR is driven largely 
by a decrease in the number of non-detections at high- 2 :, 
which we will discuss next. 


3.2. The Evolving Fraetion of Star-Forming BCGs 

In Figure |4] we show how the fraction of BCGs with 
SFR > 10 Mq yr“^ (hereafter fsr) has evolved over 
time. Within the sample presented here, this fraction 
evolves from ^20% at 0.25 < ^ < 0.65, to ^50% at 
0.65 < 2 : < 1.05, to ^90% at ^ > 1.05. This evolution 
is visible in Figure [3] as a relative lack of non-detections 
at the highest redshifts compared to the lowest redshifts. 
This result indicates that the fraction of clusters har¬ 
boring a starburst in their central, most-massive galaxy 
grows by a factor of ^3 between 2 : ^ 0.6 and 2 : ^ 1.1. 

Even more intriguing is that, at 2 : ^ 0, r\j 0 - 10% 
(jPonahne et al.l[201(l iFraser-McKelvie et al.ir2014l ). sug¬ 
gesting an even more dramatic decline since z > 1. 
We compile in Figure |4] estimates of fsF from various 
surveys, spanning 0 < 2 : < 0.5. These su rveys in¬ 
clude both optically-selected (jMcDonaldll201lL McDon¬ 
ald e t al. i n prep) and X-ray-s e lected (iDonahue et ah 
1992 . I 2 OIOI: iHoffer et al.l 120121 iFraser-McKelvie et al. 
2 OI 4 I iDonahne et al.ll20"T5l ) clusters, and are being com- 
pared here to an SZ-sel ected sample. Fo r the ACCEPT 
(iCavagnolo et al.l l2009f) and CLASH (jPonahne et al.l 
l2015f ) samples, we apply a small correction due to the 
fact that both of these cluster samples are biased to¬ 
wards cool core clusters. This bias correction assumes 
that the tru e, underlying fraction of cool core clus¬ 
ters is 30% (iHaarsma et ^ 120101: [Hudson et al.l 120101: 
iMcDonald et al.l [201301 ) and that no n-cool cores do not 
have star-forming BCGs at ^ < 0.5 (lO’Dea et al.l [20081: 

ICavagnolo et al.ll2008l: iMcDonald et aL I I2010I ). This well- 

motivated correction brings the observed value of fs f in 
line with other, more representative samples. Encourag¬ 
ingly, at 2 : ^ 0.4 there are five distinct measurements of 
fsF^ based on three different selection methods, which all 
agree that the fraction of star-forming BCGs is 20 ± 5%. 
The agreement between our measurement of fsF from 
SPT at 0.25 < ^ < 0.5 and these earlier works suggests 
that our methodology is sound. 

We show, in Figure IH a linear fit to the star-forming 
BCG fraction, fsF^ as a function of redshift. The best¬ 
fitting line passes through the origin and has a slope of 
0.57, implying a rapid evolution in fsF over the past ^9 
Gyr. Such a fit is not physically well-motivated, and does 
not capture the rapid growth at late times. However, re¬ 
gardless of the choice of parametrization, it is clear from 
FigurelUthat the fraction of clusters harboring a strongly 
star-forming BCG at 2 : ^ 1 is signifieantly (>3cr) higher 
than at 2 : ^ 0. This evolution is perhaps unsurprising, 
since galaxies in general were more star-forming at early 
times. To address this comparison, we next compare the 
BCG evolution to that observed in field galaxy and other 
cluster members over similar redshift intervals. 

3.3. Speeifie Star Formation Rates of 
BCGs at 0.25 < z < 1.2 

When comparing BCGs to other galaxies, it is nec¬ 
essary to normalize by the stellar mass of the individ¬ 
ual galaxy. For a typical BCG, with ~ 10^^ Mq, a 
SFR of 10 Mq yr“^ is negligible in terms of contribut¬ 
ing to the overall mass of the cluster, requiring 100 Gyr 
to double the stellar mass. However, this same star for¬ 
mation rate in a low-mass galaxy like M82 is enough 
to power massive outflows and modify the galaxy mor- 
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Fig. 3. — Left: Measured star formation rates (filled circles) and upper limits (arrows) for the BCGs in this sample as a function of 
redshift. We show individual measurements from each method (UV continuum, IR continuum, [Oil] emission line) where available. The 
horizontal dashed line represents our cutoff depth of SFR = 10 Mq yr“^. We highlight the most restrictive upper limit for each BCG 
with a black outline, showing that the best constraints at high- 2 ; come from deep rest-frame UV imaging, while at low- 2 ; archival data 
from WISE is able to provide robust upper limits. Right: Similar to left panel, but now showing only the most constraining upper limit 
or detection (red) for each BCG. Where multiple SFRs were measured for a given cluster (see Figure [T]), we take the average. Systematic 
uncertainties (see §2.5) are shown as vertical error bars. For UV- and [Ollj-based SFRs, the uncertainty is dominated by the extinction 
correction, while for IR-based SFRs the uncertainty is dominated by the extrapolation from observed 22iim to rest-frame 24p,m. Black 
rectangles show the average SFR in different redshift bins, where the height represents the combined statistical uncertainty in the mean 
and additional uncertainty due to non-detections. 
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Fig. 4. — Fraction of BCGs with SFR > 10 M© yr“^ {/sf) as 
a function of redshift. Black points show data from this work in 
four different redshift bins, while col ored points show dat a from 
previous works over 0 < 2 ; < 0.5 JDonahim^e^^yj^OOJ, 

” IcDonak 3 [20nl: IHofter et al.l I Fraser-McKeMeet all 


2010 


2014 


IDonahue et iDEoIS). Samples marked with an asterisk in the leg¬ 
end have had a bias correction applied, assuming that non-cool 
cores have passively-evolving BCGs and that the true underlying 
fraction of cool cores is 30% (see §3.2). This plot demonstrates 
consistency between a wide variety of samples based on X-ray, op¬ 
tical, and SZ selection. We find a steady rise in the fraction of 
star-forming BCGs, from fsF ~ 0% at 2 ; ~ 0 to fsF ~ 70% at 
2 ; ~ 1.2. This growth is depicted with a best-fit dashed line. Ver- 
tical erro r bars are derived based on the methods described by 
ICameronI (120111) for binomial populations. 


TABLE 1 

Average BCG Properties 




z 

Nei 

(M,) 

[1012 Mo] 

(SFR) 
[Mq yr-i] 

{sSFR) 

[Gyr-i] 


0.25 

- 0.45 

19 

2.09 

4.0 - 15 

0.002 - 0.007 


0.45 

- 0.65 

28 

1.97 

7.0 - 24 

0.004 - 0.012 


0.65 

- 0.85 

22 

1.47 

20 - 41 

0.013 - 0.028 


0.85 

- 1.05 

8 

1.39 

42 - 84 

0.030 - 0.061 


1.05 

- 1.25 

6 

1.54 

89 - 320 

0.058 - 0.205 


Note. — Average values for BCGs in five redshift bins, 
as shown in Figures [3] and [5] The ranges quoted on (SFR) 
and (sSFR) correspond to the combined uncertainty in the 
measurements (bootstrapping errors), the uncertainty in the 
BCG choice, and the added uncertainty associated with stack¬ 
ing non-detections. The total number of clusters is less than 
the number of BCGs due to the fact that we consider multiple 
potential BCGs for several clusters. 


pholog:y and stellar content o n short (Myr) timescales 
(iForster Schreiber et alJl2QQ3f ). Thus, if we want to com¬ 
pare BCG evolution to field galaxy evolution, we must 
consider instead the specific star formation rate, or sSFR, 
defined as sSFR = SFR/M*. 

In Figure [ 5 ] we show the sSFR for BCGs in this work, 
as well as BCGs in s imilar-mass, low- 2 : clu s ters f rom 
iHaarsma et alJ (j2QlQf ) ; iFraser-McKelvie et al.l (j2Q14f ) . In 
calculating the sSFR, we use the stellar mass of the old 
population obtained in the SED-fitting process (see §2.3). 
For comparison, we s how the average sSFR for field and 
cluster galaxies from [Alberts et al.l (j2Q14f ) . This earlier 
work showed that, at present, cluster galaxies have sup¬ 
pressed star formation compared to the field. However, 
at 2 ; ^ 1.2, galaxies in clusters are as star-forming as 































































their field counterparts, and are evolvingi: m ore rapidly. 
To match the analysis of lAlberts et al.l (j2Q14f ) , who calcu¬ 
late sSFR(z) by stacking far-IR data on galaxy positions, 
we estimate average sSFR values for our data by sepa¬ 
rately summing the total SFR and total stellar mass for 
galaxies in five redshift bins. For BCGs with upper lim¬ 
its on their SFR, we assume two limiting cases: the case 
where the SFR is equal to the upper limit, and the case 
where the SFR is zero. The net result of this stacking is 
shown in Figure [5] and Table [TJ where the height of the 
black boxes represents the combined uncertainty in the 
measurements (bootstrapping errors), the uncertainty in 
the BCG choice, and the added uncertainty associated 
with stacking non-detections. 

Figure [5] demonstrates that, at 0.5 < ^ < 1, BCGs are 
evolving similarly to the cluster member galaxies suggest¬ 
ing a common fueling mechanism. At low- 2 : (z < 0.5), the 
evolution of BCGs is less rapid than in the cluster envi¬ 
ronment, suggesting that the quenching processes acting 
on member galaxies (e.g., ram pressure stripping, stran¬ 
gulation) may not be affecting the central galaxy. This 
change in slope may be due to the gas reservoir being re¬ 
plenished in BCGs by cooling of the ICM, which is likely 
fueling star formation i n the lowest redshift BCGs (e.g., 
[McDonald et al.ll2011b[ ). We will return to this idea in 
the discussion below. At z ^ 1 there is marginal evi¬ 
dence that BCGs may be evolving more rapidly than the 
member galaxies, suggesting a preferential quenching in 
the cluster core. 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 
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Fig. 5.— Specific star formation rate (sSFR) as a function of red¬ 
shift for BCGs i n this work (red points, gray upper limits) and fo r 
nearby clusters (jHaarsma et al.l201CII : IFraser-McKelvie et al.l2014li . 
Black boxes show the combined sSFR in different redshift bins, 
computed as the sum of the individual SFRs divided by the sum 
of stellar masses - these points include the upper limits, with the 
vertical size of the boxes representing the associated uncertainty 
(see Table [T}. Dashed blue and purple curves show the evolution 
measured for t he fie ld and cluster environment, respectively, from 
lAlberts et al.l (I2014IV We show these curves down, for compari¬ 
son with the BCG evolution at 0.5 ^ z < 1. At low-z, the BCG 
evolution is slower than for typical cluster members, suggesting 
an additional source of fuel that the non-BCG cluster galaxies are 
not accessing, while at high- 2 ; ( 2 ; > 1) there appears to be a rapid 
increase in the amount of star formation in BCGs. 


3.4. Which Clusters Host Star-Forming BCGs? 



Fig. 6.— Cluster central entropy (Kp) as a function of redshif t 
for all clusters in this sample (see also IMcDonald et al.ll2013bl b 
Here, we separate clusters into two subsamples within each red¬ 
shift bin: those harboring the 20% most star-forming BCGs (blue) 
and those harboring the 80% least star-forming BCGs (red). At 
0. 0 < 2 : < 0.25 we show clusters from the ACCEPT database 
ofICavagnolo et al.l (120091). w here star for mation rates come from 
IFraser-McKelvie et al.l (|2014I1 ( 2 ; < 0.1) and ICrawford et HI (1199911 
{z > 0.1). For these low- 2 ; clusters, we show only systems with 
M500 > 3 X 10^"^ Mq and have recomputed the central entropy in 
an aperture with radius 15 kpc (co rresponding to ^0.01-0 .02R5oo), 
to prevent resolution bias (see e.g.. lPanagoulia et al.ll20l3l1 . In each 
redshift interval we show the mean Kp and the error on the mean 
for both subsamples, as well as the total population (black). While 
at 2 ; < 0.6 there is a clear separation in populations, with the 
most star-forming BCG s being found primari ly in clusters with 
low-entropy cores (e.g., ICavagnolo et al.l 1200811 . there appears to 
be no correspondence between cool cores at star-forming BCGs at 
^ > 0.6. 

At low-^, star-forming BGGs are, with very few excep¬ 
tions, foimbjiiMow^^trow, relaxed “cool core” clusters 
(e.g.. ICavagnolo et al.ll2QQ^ . This appea rs to still be the 
case for CLASH clusters out to 2 : ^ 0.4 (jPonahue et al.l 
2Q15D. and for individu al clusters at higher redshift (e.g. 
McDonald et al.l [2Q12c[1 . However, it has not yet been 
established whether the presence of a cool core at high- 2 : 
correlates with a star-forming BCG. 

In Figure ini we show the distribution of central entropy 
(Ko) as a function of redshift for clusters hosting the most 
star-forming BCGs compared to those hosting the most 
passive BCGs. Specifically, in each redshift bin we sepa¬ 
rately compute the average central entropy for the clus¬ 
ters harboring the 20% most star-forming BCGs (e.g., 
Perseus, Abell 1835, Phoenix, etc) and for those harbor¬ 
ing the 80% most passive BCGs. At 2 : < 0.3, we show 
clust ers from the ACCEPT sample (jCavagnolo et al.l 
1200911 . where the S FRs are inferred from arch ival in¬ 
frared data {z < 0.1: IFraser-McKelvie et al.ll2014D o r the 
Ha line luminosity (z > 0.1: iCrawford et al.l 119991 1. We 
have recomputed the central entropy in a 15 kpc aper¬ 
ture (^0.015R5 oo) foi* clusters at 2 : < 0.3 to match our 
more coa rse, high- 2 : apertme, a nd to avoid resolution 
bias (e.g.. lPanagoulia et al.ll201^ . We include only clus¬ 
ters with M 500 > 3 X 10^^ Mq (based on L^-M relation 
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Fig. 7. — This plot shows the fraction of relaxed (blue) and dis¬ 
turbed (red) clusters as a function of redshift, for two subsam¬ 
ples: clusters harboring the 20% most star-forming BCGs (upper 
panel) and those harboring the 80% least star-forming BCGs (lower 
panel). For low-redshift clusters (z < 0.3, two leftmost bins) re¬ 
laxedness h as been dete r mined based on the X-ray “symmetry” 
repo rted in IMan tz__eF_alJ (12Q15I), while star formation rates come 
from IFraser-McKelvie et ahl (120141) and fCrawford et al.l (119991 b For 
the three higher-redsh ift bins, relaxedness is quantified following 
INnrgaliev et al ] igoTS). This figure demonstrates that the most 
star-forming BCGs tend to be found in relaxed clusters at low-z, 
in agreement with the literature, while at high-z they are found 
predominantly in morphologically disturbed clusters. There is no 
strong trend in the the morphology of clusters hosting the 80% 
least star-forming BCGs. 


20% Most Star-Forming BCGs 



fueled by mergers to reside in clusters with disturbed X- 
ray morphology. In Figure 0 we test this scenario. For 
clusters in the SPT sample, we define X -ray morphology 
using the parameter, following iNurgaliey et al.l 

(I2Q13D . with a relaxed critereon of aphot < O-l- This 
quantity is less biased to signal-to-noise than other in¬ 
dicators, such as power ratios and centroid shift, but is 
consistent with these i n the limit of high sign al to noise, 
as is demonstrated by lNurgaliey et aLr(|2Q13F At low-z, 
we use the rece ntly-compiled list o f “symmetry” mea¬ 
surements from iMantz et al.l (|2Q15f ). using oyer lapping 
clusters to determine a common “relaxed” criterion. As 
before, for clusters at ^ < 0.1 an d 0.1 < ^ < 0.3, we use 
BCG star formation rates from iFraser-McKelyie et "ahl 
(|2014f ) and iCrawford et al.1 (|1999f ) , respectiyely, cutting 
on mass (M 500 > 3 x 10^^ Mq) in order to ensure unifor¬ 
mity in the samples. 

Figure 0 confirms that, at ^ < 0.3, the most star¬ 
forming BCGs tend to reside in relaxed, cool core clus¬ 
ters. The star formation in these BCGs is most likely 
being fueled by the cooling ICM, where cooling flows are 
most commonly found in relaxed clusters. On the con¬ 
trary, at ^ ^0.6, the most star-forming BCGs are found 
more often in clusters with disturbed X-ray morphology, 
with ^90% of the most star-forming BCGs at z > 0.6 
being found in such systems. This implies that, at early 
times, star formation in the BCG is more strongly cor¬ 
related with the dynamic state, rather than the cooling 
state, of the ICM - the inyerse to what is obseryed in the 
nearby Uniyerse. 

Figures [6] and 0 suggest that there is a transition be¬ 
tween star formation in low-z BCGs being linked to low- 
entropy, cool core clusters, to star formation in high-z 
BCGs being related more to a disturbed cluster mor¬ 
phology than to the cooling properties of the ICM. Be¬ 
low we will discuss possible interpretations of this result, 
and others presented thus far, while also addressing any 
potential systematic biases in this study. 


from iPratt et al.l[2QQ9[) . to mimic the S PT selection. We 
find, in agreement with preyious works (|Cayagnolo et al.l 
I2QQ8I ). that clusters hosting the most star-forming BCGs 
haye a typical central entropy of Kq ^ 30 keV cm^, while 
those hosting the least star-forming BCGs haye Kq ^ 100 
keV cm^. As the redshift increases, the ayerage core en¬ 
tropy in clusters hosting the most star-forming BCGs in¬ 
creases to ^40 keV cm^ (z ^ 0.4), and then to ^100 
keV cm^ (z ^ 0.7). In the two highest-redshift bins 
(z > 0.6), there is no statistical difference between the 
distribution of core entropy in clusters with star-forming 
and passiye BCGs. We find no strong eyolution in the 
core entropy of the full sample, or in the subsample of 
clusters with passiyely-eyolying BCGs. This seems to 
suggest that, while ICM cooling is likely responsible for 
proyiding the fuel for star formation in low-z clusters, a 
different mechanism is responsible for star formation in 
the high-z BCGs. 

Another mechanism for forming stars in BCGs is yia 
mergers with gas-rich galaxies. Under the assumption 
that such mergers happen most often shortly after the 
infall of a group or other massiye halo, which proyides 
an influx of new galaxies along orbits that may not be 
stable, we would expect BCGs with star formation being 


4. DISCUSSION 

4.1. Comparing X-ray and UV Morphology for 
Individual Star-Forming BCGs 

Based on Figures [6] and 0 there appears to be an eyoly- 
ing connection between star formation in the BCG and 
the host cluster morphology. In order to inyestigate the 
fuel source for this star formation (i.e., cooling, gas-rich 
mergers, etc), and determine whether this is linked to 
the dynamical state of the cluster, we require rest-frame 
UV imaging at significantly higher angular resolution. 
Fortunately, many of the high-z clusters in this sample 
haye been obseryed by the Hubble Space Telescope (HST) 
as part of yarious weak-leasing programs (Pis: Stubbs, 
High, Schrabback). These obseryations are all in the 
F606W filter, which corresponds to rest-frame i^-band at 
z ^ 0. 7 and <3000A at z > 1. These data proyide a 
detailed yiew of the young stellar populations in a sub¬ 
sample (10/14) of high-z star-forming BCGs. 

In Figures [8] and [9] we highlight the X-ray and UV mor¬ 
phology of the cluster and BCG, respectiyely, for all 10 
clusters with ayailable rest-frame i^-band (or bluer) HST 
imaging at z > 0.7. For each cluster, we also show a near- 
IR image of the core region, to demonstrate that these 
star-forming galaxies do indeed haye the highest stellar 
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Fig. 8. — Upper row: Smoothed X-ray images from Chandra of five clusters from this paper at 2 : > 0.7 with rest-frame near-UV imaging 
from HST and with both relaxed X-ray morphology and a low-entropy core. The field of view for each cluster corresponds to R 500 on 
a side. Middle row: Near-infrared image (FourStar or Spitzer) of the inner region of the cluster. The footprint of this image is overlaid 
on the X-ray image above. In all panels, the galaxy identified as the BCG is at the center of the field of view. Lower row: Rest-frame 
near-UV image of the BCG from HST. The footprint of each image is shown on the near-IR image directly above it. These images show 
that, for the most part, the relatively smooth UV emission in these galaxies is tracing the underlying old stellar populations, contrary to 
the highly-asymmetric, filamentary star formation observed in the cores of low- 2 :, relaxed galaxy clusters.. 



Fig. 9.— Similar to Figure [8] but now showing five clusters at 2 ; > 0.7 with disturbed X-ray morphologies. The BGGs in these clusters 
all have highly asymmetric UV emission, suggesting a different mode of star formation than the BGGs in relaxed clusters. 
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mass of all the galaxies in the core, consistent with being 
the BCG. We have divided these 10 clusters by X-ray 
morphology into two subsamples of relaxed {aphot <0.3 
and Ko < 95 keV cm^) and unrelaxed {aphot > 0.3 
or Ko > 95 keV cm^). The five most relaxed clusters 
(Figure [8]), in general, have relatively smooth, symmet¬ 
ric UV morphologies that resemble the underlying old 
stellar (near-IR) distribution. We caution that much of 
this near-UV emissio n may originate in the old stellar 
populations (see e.g. IHicks et a] .112(11(1 - we require ei¬ 
ther far-UV imaging or equally high angular resolution 
imaging in a near-IR band (allowing subtraction of the 
old population) to determine the morphology of the ex¬ 
cess UV emission due exclusively to young stars. Two of 
the BCGs shown in Figure [8] (SPT-CLJOOOO-5748, SPT- 
GLJ2043-5035) appear to be in the midst of major merg¬ 
ers, based on both the near-IR and near-UV imaging, de¬ 
spite the fact that these two BCGs reside in the two most 
relaxed high- 2 : clusters in our sample. None of the star¬ 
forming BGGs in this high- 2 :, relaxed cluster subsample 
show evidence for extended, asymmetric filaments of star 
formation wh ich are commonly found in analogous low- 2 ; 
systems fe.g..lO’Dea et al.ll2QlQl: [McDonald et al.ll2QlTbl : 

iDonahne et al.l[2QT^ Tremblay et al. 2015). Instead, the 

UV emission is concentrated in the BCG center for 4 
out of 5 BGGs, perhaps indicating that these systems 
are experiencing either nuclear starbursts or are AGN 
misidentified as star-forming galaxies. 

On the other hand, in the five least relaxed clus¬ 
ters (Figure [9]) the UV emission in and around the 
BCG is clumpy and filamentary. In particular, SPT- 
GLJ0102-4915 and SPT-GLJ2106-5844 exhibit clumpy 
UV emission extending 61 kpc and 57 kpc from the 
BCG center, respectively. This is comparable to 
the most extended star-forming filaments foun d in 
Abell 1795 (50 k pc: iMcDonald fc Veillen3 1200911 and 
Perse us (60 kpc; iGonselice et al.l 120011: iGanning et al.l 
|2014[ ). In all five of these BCGs, there are a minimum of 
two distinct UV peaks - in several systems, the brightest 
UV peak is offset from the BGG nucleus. Unlike the five 
BCGs in the more relaxed clusters, none of these sys¬ 
tems are consistent with being purely AGN, given their 
complex morphology. 

While incomplete, these HST data suggest a quali¬ 
tative difference in UV morphology between BCGs in 
relaxed, cool-core clusters and those in unrelaxed sys¬ 
tems. In relaxed, cool-core clusters, where low-^ surveys 
tend to find star-forming BCGs, we find (in 4/5 cases) 
smooth, centrally concentrated UV morphologies. In the 
more morphologically-disturbed systems (based on X-ray 
imaging), we find (in 5/5 cases) clumpy or filamentary 
UV emission, with several (>3) distinct emission peaks. 
Follow-up studies involving deep, high spatial resolution 
imaging at both far-UV and near-IR of a larger sample 
of high- 2 : BGGs are necessary to determine whether this 
emerging trend is merely coincidence or evidence for a 
link between central star formation and the dynamical 
state of the cluster core. 

4.2. An Evolving Fuel Supply? 

The results presented in §3 suggest that, at high- 2 :, the 
potential for star formation in the BCG is maximized 
when that BCG belongs to a dynamically-active cluster. 


This is contrary to the established wisdom, based on nu¬ 
merous studies of low-^ clusters, that the most relaxed, 
cool core cl usters tend to harbor t he most star-forming 
BCGs (e.g.. iCrawford et al]ll999l: lEdwards et al.l[2QQ7[ 

iCavagnolo et al.l[2QQ8l : IDonahne et al.ll2QIQ[) . We proceed 

with a discussion of these results under the assumption 
that the X-ray morphology traces the dynamical state of 
the cluster, or, more specifically, that morphologically- 
disturbed clusters in the X-ray are undergoing (or have 
recently undergone) a major merger. 

At early times, the cores of galaxy clusters contained 
a higher fraction of star-forming galaxies than they do 
today - this is known a s the “Butcher-Oemler Effect” 
(|Bntcher fc 0eml^ll984f ). This trend continues with in¬ 
creasing redshift, such that at ^ ^ 1.4, the star formation 
rate in the field and in the cluster environment are in¬ 
distinguishable (iBrodwin et al.ll2QI3 l : [Alberts et al.ll2QI^ 
IWagner et al.l 120151 . see also Figure [5]). Thus, we would 

expect that there is a significantly higher fraction of gas- 
rich galaxies in our high-^ subsample than in the low- 2 : 
systems, particularly in the inner cores where gas deple¬ 
tion is most effective at low- 2 :. Regardless of the impact 
angle of a cluster-cluster merger (as long as it is a bound 
system), the cores will typicall y pass through one an¬ 
other within ^2 crossings (e.g., iRicker fc Sarazinll2QQll: 
iPoole et al.ll2QQ6h . The crossing of these two cores will 
result in both an increase in the effectiveness of ram- 
pressure stripping, due to the high relative velocities of 
the cores, and an increased rate of galaxy-galaxy mergers 
and harassment. Each of these elevated processes result 
in the removal of cool gas from member galaxy halos. 
Assuming that this gas remains cool, this would lead to 
an increase in the availability of fuel for star formation 
in the cores of merging clusters. This is not the case in 
low- 2 : systems, since galaxies in the cores of low- 2 : clus¬ 
ters tend to be gas-poor (hence the low star formation 
rates), so very little gas is available to be removed from 
their halos and contributed to the core. 

With this scenario in mind, we reexamine Eigure [5l 
modeling the observed evolution in BGG star forma¬ 
tion with two distinct evolutionary components. These 
models, shown in Eigure [101 assume that the evolution 
can be described as the sum of two power laws in time 
{{sSFR) oc e^/'^). The best-fitting models have decay 
times of 4 Gyr (low-z) and 0.7 Gyr (high- 2 :), compared 
to 1.5 Gyr and 2.2 Gy r for high- 2 : cluster an d field galax¬ 
ies, respectively, from [Alberts et all (j2014f ). This model 
is meant to provide a qualitative assessment of what may 
be driving star formation in BCGs over time - there is 
no shortage of models with an equal number of free pa¬ 
rameters that would provide an equivalently good fit to 
these data. 

The models shown in Eigure [5] were chosen to illustrate 
two epochs of declining star formation. Erom 2 : ^ 1.4 to 
2 : ^ 0.8, the sSER declines more rapidly (r = 0.7 Gyr) 
than both the field (r = 2.2 Gyr) and cluster (r = 1.5 
Gyr) environment, implying that the BCG is quenched 
more rapidly than the other cluster members. This is 
consistent with an extra polation from t he fie ld to the 
cluster environment from [Alberts et al.[ (j2014l ) , indicat¬ 
ing that (unsurprisingly) quenching of star formation is 
a strong function of local galaxy density. At late time, 
from ^ ^ 0.6 to 2 : ^ 0, the sSER of BCGs is evolving 
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Fig. 10.— Here, we repeat Figure [5] now with a two-component 
fit to the BCG evolution. This fit assumes two epochs of star 
formation in BCGs: low-level star formation at present day, which 
correlates with the presence of a relaxed morphology and cool core 
in the ICM, and a rapid increase in star formation at 2 : > 0.8, which 
correlates with the presence of disturbed ICM morphology. We 
propose that this is indicative of two epochs of BCG star formation: 
merger-driven at early times, and cooling-induced at late times. 


more slowly (r = 4 Gyr) than both the field and cluster 
environment. Considering that field galaxy evolution is 
approximately passive, this implies that BCGs have an 
additional source of fuel - presumably the cooling ICM 
in relaxed systems. The non-zero slope of this evolu¬ 
tion may be due to an ever-improving balance between 
ICM cooling and AGN feedback in the cores of galaxy 
clusters, leading to a reduction in the efficiency of ICM 
cooling over time. 

In summary, there appears to be a transition in the 
source of fuel for star formation in BCGs over the past 
^10 Gyr. In nearby clusters, star formation is likely to 
be fueled by the cooling ICM, and regulated by feed¬ 
back from the central AGN. At early times, the most 
star-forming BCGs are found in dynamically unrelaxed 
clusters, suggesting that star formation may be predomi¬ 
nantly fueled by interactions with other galaxies, similar 
to the other, non-BCG members. 

4.3. Lowering the Preeipitation Threshold in Mergers 

Several recent studies have suggested that the condi¬ 
tion for thermal instability in the hot ICM can be de¬ 
scribed as tcooi/tff = 10, where tcp oi is the cooling time, 
and iff is the local free-fall time ( McConrt et al.ll2Q12l: 

iSharma et al.ll2Ql^lGasDari et al.ll2Q12l:IVoit fc Donahna 

I2Q15I: I Vo it et ^l2Q15[ ) . This threshold can be interpreted 

as the ratio of the cooling time to the mixing time. As¬ 
suming that AGN feedback is anisotropic, preventing 
local thermal instabilities requires that the heated gas 
mixes with the cooling gas on timescales shorter than 
the local cooling time. 

In order to cross this threshold, one of two conditions 
can be met. The gas can cool, reducing the local cooling 
time, tcooi- In the precipitation-driven feedback scenario 


of IVoit et al.l (j2Q15f ). this leads to rapid precipitation of 
cold clouds, which fuel AGN feedback, leading to an in¬ 
crease of tcooi- This cycle can repeat indefinitely, until a 
more energetic process drives tcooi to much larger values 
(^1 Gyr). Alternatively, if the free-fall time is increased, 
gas at a fixed density and temperature will more readily 
condense out of the hot phase. This criteria is, in princi¬ 
ple, met if the dense core of a galaxy cluster is dislodged 
from the minimum of the dark matter potential, as is the 
case during an interaction with another massive group or 
cluster. It remains unclear, however, whether increasing 
the free-fall time of the gas in this way would, in fact, 
lead to more favorable cooling. Strictly speaking, a dis¬ 
lodged cool core would have a lower value of tcooi/tff^ 
but the simulations which arrived at the threshold of 
tcooi jtff = 10 for thermal instability assumed relaxed 
clusters. In the case of a relaxed cluster, tff is a proxy 
for the convective timescale of the hot gas. However, in 
a merging cluster, bulk flows will likely drive mixing on 
faster timescales, perhaps leading to less efficient cooling. 

This represents an alternate scenario which would lead 
to a link between enhanced star formation in the BGG 
and an unrelaxed dynamical state of the cluster core 
at early times. We lean toward the merger-induced SF 
scenario to explain the uptick in star-forming BGGs at 
early times, but stress that follow-up studies, including 
spatially-resolved far-UV (e.g., HST) and far-IR (e.g., 
ALMA) imaging and optical spectroscopy, of these BCGs 
are necessary to provide additional insights into the 
mechanism for enhancing star formation in these high-z 
BCGs. 


4.4. AGN Contamination 

As we discussed in §2.6, it is possible that many of 
the BCGs that we have identified as star-forming may, 
instead, be AGN. With the inhomogeneous data set in 
hand, it is challenging to differentiate between starburst 
and AGN, or composite systems, in a uniform way for 
the full sample. Thus, it may be that both the estimated 
fraction of star-forming BCGs and the absolute star for¬ 
mation rates of these BCGs may be biased high. 

Considering only clusters at 2 ; > 0.75 (the two highest 
redshift bins in Figure |4]) we detect star formation in 9- 
10 out of 22 (41-45%) of BCGs, where the uncertainty 
represents cases where multiple BCGs were identified in 
a given cluster. This is a factor of two higher than the 
fraction of star-forming BCGs at 2 ; ^ 0.4. If we con¬ 
sider only the subsample of clusters with HST follow-up, 
and label systems with symmetric, nuclear UV emission 
as non-star forming (see Figures [8] and [9)), we find evi¬ 
dence of extended, asymmetric star formation in 4 of 10 
(40%) of BCGs. Given the reduction in sample size by 
requiring HST follow-up, the uncertainty on this frac¬ 
tion is larger, with the 95% confidence interval being 
16.7-69.2%. While this is consistent with th e observed 
star-forming fraction of ^ 18% at ^ ^ 0.4 (jMcDonaldl 
I 2 OIII : iDonahue et al.l l2015f ) , it is incons istent with the 
lower value of ^5% measured at 2 ; ^ 0 (|Donahne et al.l 
I 2 OIOI: l^aser-McKelvie et al.ll2014j) . Thus, while we can 
say with high confidence that the fraction of star-forming 
BCGs has evolved significantly from 2 : ^0 to 2 ; >0.4, 
our inability to differentiate between central AGN and 
nuclear starbursts limits our ability to say with certainty 
whether this fraction continues to grow. 
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It is important to note, however, that the fraction of 
nearby BCGs harboring AGN that are bright in the UV 
or IR continuum is low - the vast majority of cluster¬ 
centric AGN are radio galaxies which appear quiescent 
at most other wavelengths. If the high fraction of sys¬ 
tems identified here as star-forming are instead AGN, 
it implies a high accretion rat e onto the superma ssive 
black hole at early times (e.g., iRussell et al.l l2Ql^ . Of 
the hundreds of known BGGs at low-z, there are only 
a handful of systems that appear to ha rbor rapidly- 


accreting, radiatively-efhcient AGN fe.g.. iRnssell et al. 

2010 

: 10 ’Sullivan et al. 

120121: lUeda et 0111201.11: IKirk et al. 

2014 

: iRnvnolds et al.l 

2014 IWalker et al.ll2014f). Thus. 


we can say with confidence that there is a dramatic in¬ 
crease in the amount of “activity” in BGGs from ^ ^ 0 
to 2 ; ~ 1 - whether that activity refers to massive bursts 
of star formation or the rapid growth of central super- 
massive black holes remains an open question. 

5. SUMMARY 

We present multiwavelength observations and inferred 
star formation rates for 90 brightest cluster galaxies 
(BGGs) in SPT-selected galaxy clusters at 0.25 < z < 
1.25, all of which have archival X-ray data. The main 
results from this study can be summarized as follows: 

• We find a significant number of BCGs (31/90) with 
SFR >10 Mq yr“^, representing a much higher occur¬ 
rence rate than that observed in galaxy clusters at 2 : ^ 
0 (^ 1-5%; iDonahue et al.l [20T(1 iFraser-McKelvie et al.l 
l2014f ). Of these 31 BCGs, one quarter (7/31) have SFR 

> 100 Mq yr“^. 

• The fraction of clusters harboring a star-forming (SFR 

> 10 Mq yr“^) BCG is found to be ^20% at 2 : ^ 0.4, 
consistent with many earlier works based on optical and 
X-ray cluster selection. This fraction rises rapidly at 2 : > 
0.8, to a measured value of ^90% at 2 : ^ 1.1. 

• The specific star formation rate {sSFR = SFR/M^) 
of BCGs has evolved more slowly (r = 4 Gyr) in recent 
times {z < 0.6) than the overall cluster (r = 1.5 Gyr) 
and field (r = 2.2 Gyr) populations. This is most likely 
due to a replenishment of gas in the BCG via cooling of 
the ICM. At early times, the evolution was more rapid 
(r = 0.7 Gyr), with the sSFR in BCGs dropping from 
^0.1 Gyr“^ at 2 : = 1.2 to ^0.02 Gyr“^ at 2 : = 0.8. 

• At 2 : > 0.6 there is no significant correlation between 
the central entropy of the host cluster and the presence of 
star formation signatures in the BCG, contrary to what 
is observed in nearby clusters. 

• While, at 2 : ^ 0, star-forming BCGs are found in 
the centers of relaxed, cool core clusters, this trend ap¬ 
pears to reverse at high- 2 :. At 2 : > 0.6, the most star¬ 
forming BCGs in this sample are found in the cores 
of morphologically-disturbed clusters (based on X-ray 
asymmetry). 

• Excluding the Phoenix cluster, the most strongly star¬ 
forming systems in this sample have SFRs of order ~100- 
300 Mq yr“^. Based on rest-frame near-UV follow-up of 


a subsample of high- 2 : systems with HST, we find that 
this star formation can be extended on scales of ^50-60 
kpc. 

The observation that an enhancement in BCG star for¬ 
mation correlates with the dynamical state of the cluster 
at high-^ suggests that star formation may have been 
fueled by interactions with gas-rich satellites at early 
times. Further studies, utilizing deep, high angular reso¬ 
lution far-IR and far-UV imaging and integral field spec¬ 
troscopy of a larger sample of BCGs will help determine 
if such mergers are, indeed, the dominant source of star 
formation in BCGs at early times. 
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Fig. A.l.— Broadband spectral energy distribution for 90 BCGs in this sample. Wavelengths here are in the observed frame. We overplot 
the best-fitting models from 100 fits, while varying unknown quantities such as the intrinsic reddening {E(B — V) = 0.3 zb 0.1) the mid-IR 
slope (a = 2.0 ± 0.5) and the formation epoch of the old population (zf =2 — 5). The model consists of an old stellar population, formed 
at z = Zf (fLeitherer et al.lfl999lb a young stellar population which has been continuously forming stars for 30 Myr (|Leitherer et al.iri999lb 
and a dust component parametrized by a powerlaw in the mid-lR, following ICasevI (|2012l i. 
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BCG Star Formation Rates and Stellar Masses 
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2 1 7+0-16 
^•4'-0.14 

- 

<5.5 

-15+4.4 

40-4.4 

SPT-CLJ0411-4819 

0.424 

62.8154 

-48.3175 


<12 

- 

27tf° 

SPT-CLJ0417-4748 

0.581 

64.3461 

-47.8132 

Q qqTO.23 

^•'9'9-0.63 

<93 

<3.8 

<28 

SPT-CLJ0426-5455 

0.630 

66.5171 

-54.9253 

2 89+0-07 
Z.OU_i 20 

<4.7 

- 

51I22 

SPT-CLJ0438-5419 

0.421 

69.5734 

-54.3224 

3 98+®'^^ 
'9.90-0 20 

<31 

<23 

<24 

SPT-CLJ0441-4855 

0.790 

70.4497 

-48.9233 

1 38+0-4^ 
i.oo_0 23 

<13 

- 

<61 

SPT-CLJ0446-5849 

1.186 

71.5157 

-58.8304 

1 46+0-09 
4-^'9-0.59 

2001^60 

- 

330±?^o 

SPT-CLJ0449-4901 

0.792 

72.2819 

-49.0214 

1 36+0-04 

4-'9'9-0.19 

<10. 

- 

<84 

SPT-CLJ0449-4901 

0.792 

72.2669 

-49.0276 

1 oi +0.14 

4-^4_o 22 


<27 

<81 

SPT-CLJ0456-5116 

0.562 

74.1171 

-51.2764 

1 48+0-13 
4-^°-0.34 

<79 

<6.3 

<32 

SPT-CLJ0509-5342 

0.461 

77.3393 

-53.7035 

1 54+0-44 
4-0^-0.16 

<32 

33120 

<16 

SPT-CLJ0516-5430 

0.295 

79.1556 

-54.5004 

2 90+°'^® 
^•4''9_o.17 

<13 

<27 

9 5+1-3 
^•^-0.82 

SPT-CLJ0528-5300 

0.768 

82.0222 

-52.9981 

1 59+0-06 
4-0^-0.13 

<10. 

<15 

<34 

SPT-CLJ0533-5005 

0.881 

83.4033 

-50.0958 

n 75+0-03 

O'- ' O_0.08 

<16 

50±?9 

<130 

SPT-CLJ0542-4100 

0.642 

85.7085 

-41.0001 

1 86+°-^^ 
i.ou_0 16 

<1.4 

- 

<28 

SPT-CLJ0546-5345 

1.066 

86.6573 

-53.7588 

1 70+0-09 

4-''9_o.30 

<8.0 

- 

iioIm 

SPT-CLJ0551-5709 

0.423 

87.8931 

-57.1451 

0 69+®'®® 

0'-'9^-0.02 

<37 

- 

<17 

SPT-CLJ0559-5249 

0.609 

89.9301 

-52.8242 

3 33+0.10 

4-4J--0.ll 

<11 

<16 

<19 

SPT-CLJ0655-5234 

0.470 

103.9760 

-52.5674 

0 36+0-00 

0'-'9'9_o.02 

<3.1 
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TABLE A.l 

BCG Star Formation Rates and Stellar Masses (continued) 


Cluster 

^cluster 

olbcg 

Sbcg 

M* 

SFRc/y 

SFR[o/ji 
[M© yr-l] 

SFR24/j,m 



[“] 

[“] 

[1012 Mq] 

[M© yr-i] 

[M© yr-l] 

SPT-CLJ2011-5725 

0.279 

302.8620 

-57.4196 

0 +0.02 
' ^-0.10 

- 

- 

- 

SPT-CLJ2031-4037 

0.342 

307.9720 

-40.6252 

1 qq+0.02 

<51 

- 

<7.5 

SPT-CLJ2034-5936 

0.919 

308.5390 

-59.6042 

1 Q7+0-18 
-0.30 

<4.5 

- 

<99 

SPT-CLJ2035-5251 

0.528 

308.7950 

-52.8564 

n fTo+0.05 
' ^-0.01 

<2.8 

<2.4 

<41 

SPT-CLJ2035-5251 

0.528 

308.7930 

-52.8539 


120“^^® 

- 

q-| +29 

0 J--I 6 

SPT-CLJ2043-5035 

0.723 

310.8230 

-50.5923 

0.88l°|6 

210“^^'^® 


<89 

SPT-CLJ2106-5844 

1.132 

316.5190 

-58.7411 

0.83l°J« 


- 

<260 

SPT-CLJ2135-5726 

0.427 

323.9060 

-57.4418 

-0.03 

<13 

<8.5 

<13 

SPT-CLJ2145-5644 

0.480 

326.4660 

-56.7482 

9 qn+0.05 

<3.6 

<26 

- 

SPT-CLJ2146-4633 

0.933 

326.6470 

-46.5505 

i.yO-o 25 

<9.6 

<34 

- 

SPT-CLJ2148-6116 

0.571 

327.1780 

-61.2795 

-*■•^^-0.05 

<7.5 

- 

<39 

SPT-CLJ2218-4519 

0.650 

334.7470 

-45.3145 

1 1 9+0.03 
-^•-^^-0.09 

4 4 + 5.0 
^•^- 2.4 

- 

<55 

SPT-CLJ2218-4519 

0.650 

334.7500 

-45.3162 

i.UD-o 27 

< 1.0 

- 

- 

SPT-CLJ2222-4834 

0.652 

335.7110 

-48.5764 

1 47+0-05 
-0.08 

<3.1 

- 

<52 

SPT-CLJ2232-5959 

0.594 

338.1410 

-59.9980 

9 9/y+0.38 
-0.26 

<110 

<18 

<33 

SPT-CLJ2233-5339 

0.480 

338.3150 

-53.6526 

9 [;:o+0.27 

Z.00-0,20 

<86 

<12 

<13 

SPT-CLJ2236-4555 

1.162 

339.2140 

-45.9295 

o.zy -0 97 

30+00 

- 

- 

SPT-CLJ2245-6206 

0.580 

341.2590 

-62.1272 

/.uy-o 86 

<5.2 

- 

- 

SPT-CLJ2258-4044 

0.826 

344.7010 

-40.7418 

^•'^^-0.07 

10 

- 

- 

SPT-CLJ2259-6057 

0.750 

344.7540 

-60.9595 

9 -| c^+0.44 
^•-'-^-0.25 

<11 

- 

<84 

SPT-CLJ2301-4023 

0.730 

345.4700 

-40.3868 

1 04+°-°^ 

36+53 

00-19 

- 

- 

SPT-CLJ2306-6505 

0.530 

346.7230 

-65.0882 

1 34+0.06 

<14 

< 2.6 

<39 

SPT-CLJ2325-4111 

0.358 

351.2990 

-41.2037 

'^•'-'^-0.13 

<9.8 

<11 

<6.7 

SPT-CLJ2325-4111 

0.358 

351.3000 

-41.1991 

i.yU-0,03 

<14 

<5.0 

<8.4 

SPT-CLJ2331-5051 

0.576 

352.9630 

-50.8650 

1 [;:q+0.06 

i.OO-0,08 

9q+36 

^0_i4 

- 

<46 

SPT-CLJ2335-4544 

0.547 

353.7850 

-45.7391 

1 40+0-22 

41 

^-'■-27 

<6.7 

- 

SPT-CLJ2337-5942 

0.775 

354.3550 

-59.7058 

0.13l°;°ii 

1 00 “^^'^ 
l.UU_o 66 

- 

- 

SPT-CLJ2337-5942 

0.775 

354.3650 

-59.7013 

2 . 081 “;°^ 

<4.9 

<40 

<76 

SPT-CLJ2341-5119 

1.003 

355.3010 

-51.3291 

1 13 + 0.10 

<2.5 

- 

<170 

SPT-CLJ2342-5411 

1.075 

355.6910 

-54.1847 

u.yu_o 22 

00+42 

00-19 

<17 

<240 

SPT-CLJ2344-4243 

0.596 

356.1830 

-42.7201 

o.yo _2 43 

2000 “^^'^°° 

^OUU_3^200 

17001“°° 

2000 l®|o 

SPT-CLJ2345-6405 

0.937 

356.2510 

-64.0927 


<71 

- 

250<f 

SPT-CLJ2352-4657 

0.734 

358.0680 

-46.9602 

-‘-•'^^-0.14 

<9.5 

- 

- 

SPT-CLJ2359-5009 

0.775 

359.9280 

-50.1672 

1 93+0.05 

0 9+13 
^•^- 6.0 

<21 

<84 


Note. — Positions and star formation rates of BCGs used in this study. UV- and [O ii]-derived SFRs 
are corrected f or int rinsic extinction assuming E{B — V) = 0.3 ± 0.1 and a gray extinction curve from 
ICalzetti et a D <20001) . Systematic uncertainties, which are dominated by the extinction correction in the 
blue bands and extrapolation to 24p,m in the IR band, are quoted for all detections. These uncertainties are 
discussed in detail in §2.5. 









